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1. Introduction 



Holographic techniques based on the AdS/CFT correspondence [1, 2] have become a fruitful arena 
in which to address questions concerning properties of a strongly coupled non-abelian plasma (see 
[3] and references therein). The finite temperature scenario [4] shares many more properties with 
what is expected from thermal QCD than its zero temperature counterpart. Within this context, 
holographic spectral functions [5] have been the subject of much attention in the last year [6, 7, 8, 10]. 
They carry information about transport coefficients such as the conductivity and viscosity of the 
plasma [11] as well as about particle production rates [12]. This is very appealing, since photon 
and dilepton production are amongst the most interesting signatures of a quark-gluon plasma. In 
fact these two phenomena have been studied in a number of recent papers starting with [13] where 
the geometry is that of a stack of black D3-branes, corresponding to pure A/" = 4 SYM at finite 
temperature with a weakly gauged U{1) r subgroup of SU (4)r, under which the Af = 4 fermions and 
scalar fields are charged. It was shown that in order to calculate the photon and dilepton emission 
rates, at small electromagnetic coupling, the dependence on the U{1) gauge field is subleading in the 
expression for current-current correlators and, therefore, the calculation can be performed solely in 
reference to the non-abelian gauge fields, with no need for a dynamical photon. This simplifies the 
problem greatly as the full gravity dual of A/" = 4 SYM coupled to weakly gauged electromagnetism 
is not known. 

A step towards a more realistic model of QCD required the presence of charged fundamental 
matter. This technology is now well under control in the quenched approximation, and involves the 
introduction of Nf <C Nc probe branes in the dual gravitational background [14]. The study of this 
framework has revealed an interesting phase structure, in which the adjoint and fundamental matter 
deconfine at different temperatures [15, 16, 17]. Photoemission in such setups has been recently 
studied in [7] in the D3/D7 and the D4/D6 systems, in [18] in the context of the Sakai-Sugimoto 
model [19], and in [20] for the AdS/QCD background. 

In order to unravel the phase diagram of QCD-like theories it is mandatory to go beyond the vanish- 
ing chemical potential limit, fi = [21, 22]. Recently, the study of the strongly coupled deconfined 
plasma from AdS/CFT with nonzero baryon or isospin chemical potential has become more than 
an academic exercise [8, 9, 17, 23, 24, 25, 26, 27, 28, 9]. For instance, at low temperatures and high 
densities one expects to find interesting new phases, like the color-fiavor locked superconductor, 
that could occur in the interior of very dense neutron stars. At large Nc a number of new phases 
have been proposed in [29, 30] which may have implications for real QCD. Concerning heavy ion 
experiments, the phenomenological appeal of this extension is less clear. While in the fireball the 
residual baryon density at the core is tiny, in experiments at SPS, fits are nonetheless consistent 
with values of the chemical potential of fio ~ AOOMeV [31, 32]. Whether this will be enough to 
locate the critical point between a crossover and first order phase transition is still unclear. There- 
fore one ought to address issues like photon and dilepton emission as a means of gaining insight 
into different regions of the phase diagram. 

The present work is an extension of [7] to include both chemical potential and finite spatial mo- 
mentum. Note also that in [8, 20, 33] the authors made some progress in this direction, although 
for simplicity they studied only a subsector of the possible spectral functions. 
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Next we review the strategy of the computation, as discussed in [13]. The basic object to compute 
is the spectral function, related to the retarded two-point function G^^, as follows 

X,Ak) = -2lmG^,. (1.1) 

At zero temperature, the form of the retarded correlator is dictated by Lorentz invariance and gauge 
symmetry 

Gj, = P,Am{k') , (1.2) 

where P^^, = ri^y — k^ki,/k'^ is the transverse projector in Minkowski space and rj^y = diag(— 1, 1,1,1). 
In A/" = 4 supersymmetric Yang Mills theory, the form of n(/c) is dictated by scale invariance up to 
a constant that can be computed and yields, at large Nc, [13] 



Xl:'{k) = P,.(A:)^|P|e(-P)sgn(A:°) . (1.3) 



At nonzero temperature T ^ only rotational invariance remains unbroken, and the tensor de- 
composition of the retarded correlator now defines two polarization tensors Il-^{k) and nll(A;) as 
follows 

G^k) = P^uik) n\k) + Pl{k) U\k) , (1.4) 
with the transverse and longitudinal projectors defined as follows: Pi^ik) = 0,Pf^{k) = 0,P^^{k) = 
6ij - kikj/k"^, Pjlu{k) + P^^{k) = P^^{k) = r]^^ - k^k^/k"^. For example if /c^ = (w, g, 0, 0) we have 

P± _ P± _ 1 . pii _ . p\\ _ -^'^ . pll _ n 5) 

Observables, like the conductivity and particle production rates, are encoded in the trace of the 
spectral function, which can decomposed as 

X^'^ik) = -4Imn^(fc) - 2Imnll(A;) = x^{k) + -x"(A;) . (1.6) 

For lightlike momenta, /c^ = 0, we see that Pjlu diverges. Hence, on the light cone, we have that 

lim nil (A;) =0. (1.7) 

Otherwise we would have a divergence of G^. Therefore, x^i^ik) at lightlike momentum, is controlled 
by the transverse polarization n-*- which is enough to compute the emission rate for real photons 
[12] 

^1,3 2 

^r,= -— — n^(A;)x%(A;) . (1.8) 

Here, the Boltzmann factor, nB{k) = (e^"/^— 1)~^, receives no contribution from the baryon chemical 
potential because it refers to the bath of thermal photons which have zero baryon number. When 
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dealing with timelike momenta, rotational invariance relates these two polarizations for vanishing 
three-momentum q where k'^ = (cu, q), 




lim n (cij, q) = lim n"(ct;, q) . 



(1.9) 



For nonzero q both polarizations contribute to X^^l{k) ^i^nd encode the production of virtual photons 
which eventually decay into dilepton pairs, ZZ, of momentum + kj = k^. If the lepton Z has mass 
mi and charge ei, the differential dilepton emission rate per unit four- volume is given by 



As we see, the relevant part of the computation resides in the retarded correlator. At strong 
coupling, for the time being, only holographic techniques are easy to implement, albeit in A/" = 4 
supersymmetric Yang Mills at large N^. 

In [34] the quasinormal spectrum for strongly coupled finite temperature A/" = 4 SYM was calculated 
from a holographic perspective. The 'melting meson' scenario defines a spectrum of states where the 
poles in the correlator are found in the complex plane and therefore correspond to modes of finite 
lifetimes. For this reason, the interpretation of the spectral functions which we will be calculating 
in the following will be clear in terms of the set of quasinormal mesons in the plasma. 
The outline of the paper is as follows. In section 2 we outline the brane construction of flavor 
probes in the presence of a finite baryon number density and discuss the regions of interest in 
the parameter space in terms of the effective horizon area on the probe. In section 3 we discuss 
the calculation of the transverse gauge fluctuations on the probe brane and the coupled system 
involving the longitudinal modes and the scalar perturbations. We show how the spectral function 
can be calculated in all three sectors. In section 4 we give results in various regions of parameter 
space for the spectral function, the photoproduction rates and the limiting velocity of the mesons, 
as calculated from the peaks in the spectral functions. We illustrate how the presence of a finite 
baryon number and finite spatial momentum affect all these results. In section 5 we discuss how 
the conductivity can be calculated using the microscopic approach and show that we get the same 
answer as that calculated in the macroscopic regime in [23]. In section 6 we discuss these results 
and provide new directions for future research. In the appendices we provide more details of the 
calculations given in the bulk of the paper as well as providing some exact, analytic results in 
various regions of momentum space. 

2. Holographic setup 

In the framework of the AdS/CFT correspondence, the retarded correlator G^{k) can be obtained 
from the perturbations of a U{1) gauge field dual to the electromagnetic current on the boundary. 
The relevant holographic description is provided by an AdS geometry with a non-extremal horizon 
and embedded probe branes. The baryonic U{1) symmetry is the abelian center of the natural 
U{Nf) global symmetry present on a stack of Nf coincident D-branes. For the case of interest here. 



dT 



II — 



d^k e'^ef 
(27r)4 67r|P|5/2 



e(A;°)e(-A;2 - Am^i)^ -k^ - Amf{-k'^ + 2ml)nB{k)x^ ^{k) . (1.10) 
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namely Dp/Dq configurations, the dynamics of this gauge field is fully encoded in the action for 
the probe Dq brane: 

S=-NfTD^ [ (P+^x e-'^^- det(^ + 2-Ka'F) + WZ . (2.1) 

JDq 

The second term on the r.h.s. stands for the Wess-Zumino term which will not make any contribution 
to the equations of motion for the background and the fluctuations (see Appendix D). T^q = 
Ij {{2nlsygsls) is the Dq-brane tension, Qs is the string coupling constant and g^y is the puUback 
metric induced by the relevant background. As for the background, we will be concerned with the 
near horizon limit of a stack of non-extremal Dp-branes. The general form for any p is given by 

ds"" = H~'/\-fdt^ + dx^) + H^/^ + p^rfl^Lp) , 

e* = i/^ ; C01...P = H-^ , (2.2) 
where x = {x^, ...,x^) and 

^(P)=(^)'" ; f(p) = ^-{jj~'- (2-3) 

The probe Dq-branes wrap an n— sphere in the directions transverse to the Dp-branes, so it is 
convenient to write the metric on 5*^^^ in adapted coordinates, 

dnl_j, = de^ + sin^ e dnl + cos^ e dn^j_p_^ . (2.4) 

Setting ip = cos 6 the classical Dq-brane embedding may be specified by a dependence ip = ipip). 
On the probe brane, a U{1) gauge field can be switched on that will also depend only on the radial 
coordinate, A^{p). We shall make use of the dimensionless radial coordinate u, related to p by 

in terms of which f{u) = l — u^ and the horizon lies at m = 1. The Hawking temperature is given 
by 

7-p fpo — 



2.1 The D3/D7 system 

From this point on we will specialize to the case of D7-brane probes in a black D3-brane geometry. 
The D3/D7 intersection is summarized in the following array 

0123456789 
D3: X X X X (2.7) 

D7 : xxxxxxxx 
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and the bulk metric reads 



^,2 ^ (![^(_Jrf^2 ^ ^^2) + + ^ (2.8) 



where 

= Ang.NJt . (2.9) 
Specifiying the D7-brane embedding through ip = tp{u) the induced metric takes the form 

o^Sd7 = — — (-M +dx ) + 4^2^(1 _ ^2) + ^ ~ ^ )^^3 ' (2-10) 

with the D7-brane wrapping an G S^. The generahzation of the previous setup for finite baryon 
density was investigated in [17]. The relevant bulk degree of freedom dual to the baryon chemical 
potential is the Aq component of a U{1) gauge field on the worldvolume of the D7-brane. The 
background profiles for ip{u) and Ao{u) are obtained by solving the Euler-Lagrange equations of 
the Born-Infeld lagrangian 

C = -NjTd.^J- det{g + 2Tia'F) . (2.11) 

The gauge field Ao{u) obeys a conservation equation owing to the fact that it enters the action 
purely through its derivatives, 

where ijj{u) = a/1 — il){uY = ^va.9{u). Asymptotically in the UV region m — > we will show that 
lim^^oV^ = (see eq. (2.19) below) and therefore this equation reduces to dl^Aoiu) = which has 
the solution 

Ao{u) = jj, - au + ... . (2.13) 

By means of the holographic dictionary [23] fi is proportional to the chemical potential for the 
baryon number density, and a is proportional to the baryon number density itself. Equation (2.12) 
implies the existence of a constant of motion, D, which we normalize as follows 

- , (2.14) 



Evaluating this constant of motion at u = implies that 



4a' 

D = ^2^- (2-15) 



In terms of D the field Aq can be expressed as 



J^,{u) = -^^===- (2.16) 
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Following the discussion in [17] we can express the chemical potential as 



/x=— ^/ ^ -du, (2.17) 



and we see that in the limit of vanishing baryon density D — we obtain vanishing chemical 
potential (note that there is a region of the phase diagram for which this does not hold for sufficiently 
large quark mass [17, 27]). 

The equation for '4){u) for a generic Dp/ Dq intersection is given in eq. (A. 12), and, specializing to 
the D3/D7 case, gives 

^« — 1 ^ = • (2-18) 

Close to the boundary this equation reads du {^ip' /u) = —3ip/u^ and its solution behaves as 

^(n)~^nV2 + ^n^/^ + ..., (2.19) 

independent of the baryon density and where the constants m and c parametrize respectively the 
quark mass and something we loosely refer to as quark condensate [35, 36, 37, 38, 39], 



(O) = -IVXNfN.T'c, (2.20) 
o 

with A = Qym^c = '^T^QsNci the 't Hooft coupling. The operator (9 is a supersymmetric version of 
the quark bilinear 

O = + <l>tx$ + Mg$^<l> , (2.21) 
with X one of the adjoint scalars. A precise definition can be found in [17]. 

Equation (2.18) is a non-linear differential equation which cannot be solved analytically and there- 
fore its integration is performed numerically [15, 41]. For D = the stable embeddings il){u) fall 
into two categories, labelled "Minkowski" and "black hole". Additionally one can find metastable 
configurations corresponding to a supercooled phase. For Minkowski embeddings the D7-brane 
never enters the black hole and we have a maximum value Umax < 1 the probe branes can reach. In 
this case one has stable bound states identified with mesons whose spectrum manifests a mass gap 
and discretization [42]. This is seen in the spectral function as an infinite sum of 5-functions. For 
black hole embeddings, the D7-brane intersects the horizon, hence Umax = 1- The branes develop 
an induced horizon and all the meson resonances become unstable [34, 40]. 
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This situations changes drastically at finite baryon density. One of the main discoveries in [17] was 
that for any value oi D > there are only black hole embeddings. In fact a regular series expansion 
around the horizon m ~ 1 that solves (2.18) takes the form 

i■(-)=i■«-l ^,^°_%^%. (^--) + 0(l-uf (2.22) 

and as we see, depends solely on the limiting value of the embedding profile at the horizon i/jq = 
ip{u = 1). Note that this solution automatically satisfies the orthogonality condition of [17] due 
to the vanishing Jacobian on the horizon. This expansion can be used to specify the boundary 
conditions in order to numerically integrate ip{u) out towards the boundary, from where we can 
read off m and c from eq. (2.19). Plotting m{ipQ) in figure 1 we observe indeed that moving the 
value of tpo in the range [0,1) covers the full range of m G [0,oo). The 3-area of the induced 

p cos 6 



m 




p sin 6 



1 2 3 4 5 

Figure 1: On the left-hand side we see the dependence of m with ipQ for values of D = 20,5 and 1 (from 
left to right). As soon as D ^ 0, arbitrarily high values of m have a black hole embedding with ■0o S [0 , 1). 
On the right-hand side we plot the embeddings ^{u) for different values of ipQ ranging from 0.1 to 0.99 and 
D = l. 

horizon (per unit 3-dimensional Minkowski space volume) is controlled by i/jq 

Ah = 2Tr\TTTLy{l-ijlf/\ (2.23) 

In figure 2 we have plotted the curves of equal induced area on the D7-brane. We expect this 
quantity to govern the rough shape of the peaks of the spectral function with larger widths for 
larger induced horizons. Close to D = we see a qualitative change in behaviour below merit ~ 1-3. 
This is roughly the critical mass below (above) which we have black hole (Minkowski) embeddings 
a.t D = 0. From the form of the iso-curves, we expect a maximum effect of D on physical quantities 
for m in the range ~ (1.3,4). For m outside this range the change in induced horizon area is small 
for a given increment in D. 
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D 




1.3 2 4 8 10 12 



Figure 2: The left plot shows the normalized horizon area as a function of the quark mass and the baryon 
number density. The right hand graph is a contour plot showing lines of equal induced horizon area, labelled 
by the corresponding value of ipQ. Hence we have smaller induced horizons for larger values ofipQ. 

3. Fluctuations 

We will consider perturbations of the world-volume fields that may depend only upon the RG 
flow coordinate u and the Minkowski coordinates x^,x^. In other words, we shall consider only 
fluctuating fields that are independent of the internal coordinates wrapping the S^. 

A,iu, x) ^ A,{u) + e e-*(-^"-^^'U^(M) . (3.1) 

With this we can expand the DBI lagrangian in powers of e 

C = Co + eCi + e^C2 + ■■■■ (3.2) 

Upon imposing the equations of motion for the background fields, Ci vanishes and the linearized 
equations for the perturbations can be derived from the quadratic piece. The equations of motion 
can be found in appendix A. For q ^ the fields Ao,Ai and \E' must solve a coupled system of 
differential equations and cannot be set to zero independently. This coupling only happens for q 
and vanishes in the limit of zero D. On the other hand, the transverse excitations A± = ^2,3 do 
decouple. We will analyze the two sets of fields independently in the following section. 
It is important to note our definition of the various degrees of freedom. Although the scalar and 
the longitudinal mode of the gauge field do mix via their equations of motion, we will show that 
the spectral functions decouple in the UV (that is there is a vanishing two-point function between 
the operators sourced by the scalar and the longitudinal vector field). This means that although 
there is IR operator mixing, there is no ambiguity in the definition of the vector and scalar modes 
on the UV boundary. 



- 9 - 



3.1 Longitudinal and Scalar excitations 



This set of degrees of freedom is rather cumbersome to analyse because their equations of motion 
are coupled at first order. We have the freedom to make a gauge choice which leaves us with a 
constraint, plus three linearly independent equations relating Aq, Ai and \& (see eqs. (A. 15 - A. 17)). 
By writing linear combinations of the constraint with the equations of motion for Aq, Ai and \1/ we 
are able to find a single propagating gauge invariant combination given by the longitudinal electric 
field component: 



This field is still coupled to the scalar excitation \l/ and thus the equations of motion for these two 
degrees of freedom must be solved simultaneously (see eqs. (A. 23, A. 24)). In order to solve the 
coupled equations of motion we must investigate the boundary conditions for the coupled system. 
For retarded Green's functions, we select incoming wave boundary conditions on the black hole 
horizon. The analysis of the Frobenius expansion is given in appendix B. The selection of incoming 
wave boundary conditions is possible and from it we can numerically integrate the coupled system 
of equations simultaneously to the boundary to obtain the Fourier bulk modes "^kiu) and Zk{u) 
with k = (ci;,g,0,0). At this stage we must write the boundary action in terms of the \l/fc(M) and 
Zk{u) degrees of freedom. This is given by an expression of the form: 



Sb[Z,^-u] = -M / -—{FzzZ'^Z^k + Fz^{Z'^^,k + %Z^k) + F^^%^.k + ...) , (3.4) 



where the ellipses stand for non-derivative terms that will not contribute to the imaginary part of 
the spectral functions and 



The coefficients, F/j, are functions of the embedding solution, D, on, q and u which are best given 
in terms of the following functions 



Z = qAo + (^Ai . 



(3.3) 





(3.5) 




(3.6) 




(3.7) 



We find, in terms of the usual dimensionless ratios t» = uj/2'kT and q = g/27rT, 
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F (h ^ _ rfg^h 1 1 

' ^ q2 + 0,2^ \^ 27ra' / 27rT uo^ - \ 27ra' / 

/V;i4/l3/2(q2(J^8 _ 4/)2j2^/2^5) ^ ^^2^2^) /7^TL2^ ' 



+ ... 



F^j,q,{k, u) 



U-.0 1 fnTL^Y 



u \ 2'Ka' 



+ .... (3.8) 



From these expressions it is clear that the couphng between operators sourced by Z and \1/ in the 
UV occurs only for non-zero q and D. 

One can show that, on shell, the imaginary part of the boundary action (3.4) is independent of u. 
Despite the fact that this is of little help in computing the individual components of the retarded 
two-point function Gfj{k), we have used this as a quality check of our numerical integration. In 
fact, obtaining a constant value up to one part in 10^ supports our confidence in the accuracy of 
the solutions. 

In order to calculate the spectral function, we must take derivatives with respect to the boundary 
values of the fields which act as sources for the operators of interest. However, with the degree 
of freedom in its current form, this is not the source for the quark-bilinear operator which we are 
interested in. By looking at the solution of the field on the boundary (eq. (2.19)) we see that the 
field with the appropriate scaling to source a qq operator (and its gaugino and scalar counterparts) 
in the UV is \E'(n) = n~2\I>(ti). Changing variables to ^ we can calculate the retarded correlator in 
the Z sector which is given by Gfj{k), I,J = Z, \E^: 



GfAk) = -AT 



v 




Fz^{k,u)u 



2uK 




(3.9) 



where the ellipses stand for real contact terms that will not contribute to the spectral function. 
The boundary limit m — > corresponds to a definition of the dual theory in the UV. We see that 
in this limit the off-diagonal terms vanish and we obtain a diagonal matrix. Conversely, setting 
u = e corresponds to a definition of the dual QFT as a Wilsonian effective field theory defined at 
a finite energy scale. In this case the presence of D will induce in the RG flow a mixing of the 
corresponding quantum operators. 
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3.2 Transverse excitations 

Transverse excitations are decoupled and tlierefore simpler to deal with in the presence of D. Again, 
computing the retarded Green function involves the evaluation of the boundary action on relevant 
solutions for the bulk perturbations. The equation of motion for these perturbations can be found 
in equation (A. 14). Close to the horizon the incoming transverse mode has the form: 

A^,k{u)r^{l-u)-''^ak{u), (3.10) 

with k = {u, q, 0, 0) and ak{u) an analytic function at n = 1. The solution for A±^k{u) will depend 
parametrically on m and D (both explicitly and through tpolm, D)), uj and q. These solutions are 
to be inserted into the boundary action which is given by 

Sb[A^M = -^fJ -0^,F^{u)A'^^,iu)A^,^k{u), (3.11) 

where F±[k, u) = ftjj^Vh, with h{u) as in (3.7) of the previous section. An analysis of the differential 
equation shows that close to the boundary the exponents of the Frobenius expansion are A = 0, 1 
as usual for a bulk gauge field. In this case, the usual prescription for calculating Minkowskian 
2-point functions is implemented [43] 



Gf{k) 



-2Af\imF^{k,u) 



A±^k{u) 



(3.12) 



In section 4 we plot the results of this analysis. When comparison is possible, we find agreement 
with the results of [8, 33]. 



4. Graphical results 

In this section we give an account of the results obtained by numerical integration of the equations 
of motion. In all cases we calculate the spectral function by using the definitions given in equations 
(1.1), (3.9) and (3.12). The equations of motion are defined in terms of the parameters q, t», D 
and implicitly through the mass m. In the following we will explore what we believe to be the most 
enlightening areas of this large parameter space. In various limits, the numerical results can be seen 
to coincide with the analytical behaviour, derived in appendix C. 

4.1 Vector modes 

4.1.1 Lightlike momenta q = t» 

For the study of spectral functions at lightlike momenta, the longitudinal excitations vanish and we 
must simply solve for the transverse degrees of freedom. Spectral functions at lightlike momenta 
have been calculated for D = in [7] as they are relevant for the computation of the emission of 
thermal photons from the plasma. For lightlike momenta the longitudinal contribution to the 
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D=0 D=0.1 D=0.5 




Figure 3: In this figure we plot the quotient of the spectral function over the frecuency x^^(t»)/i» for 
lightlike momenta to = q and different values of the constituent quark mass. Above D ~ 0.3 all curves 
become monotonically decreasing with to. 



m=0 in=0.9 in=1.3 

dr dr dr 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



Figure 4: Illustration of the increase in the photoproduction rate with D for masses below merit- 



spectral function vanishes. In figure 3 we have plotted^ the spectral function for different masses 
and baryon density. The left most plot coincides with the results in [7] and the others show the 
influence of D on the set of curves labelled by different masses. 

As a general rule we find an enhancement of the spectral function for low values of the momentum. 
The top curve corresponds to massless quarks, and its behaviour at the origin is related to the 
electrical conductivity. We observe an increase with D that will be investigated below. In fact, 
for sufficiently high values of D{>^ 0.3) all the curves for different masses become monotonically 
decreasing functions of u. A peculiar feature at D = is the crossing of the curve m = 1.3 with 
curves of slightly lower mass. This, translated into emission rates suggests that plasmas made of 
heavier quarks would shine brighter beyond t» ~ 1.4. We see that this feature disappears already 
for low values of the baryon number and becomes a single valued function of u and m. The 



In all the plots wc normalize the spectral function in units of ^ — and the scalar mode in units of 
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enhancement of the spectral function, translated into photopro duct ion rates exhibits a considerable 
increase as plotted in figure 4. The white area below the coloured curves corresponds to the D = 
value. The enhancement is extremely large for masses around rricrit (the mass of the quark at the 
phase transition point in the absence of baryon number ~ 1.3). This corresponds to an effective 
increase in brightness of the plasma of more than an order of magnitude! 

4.1.2 Large n? behaviour 




Figure 5: Asymptotic behaviour of o function of the mass. We see that for large masses the 

coefficients of the fit stabilize around a value of j ~ 1. Note that numerical errors are large here due to 
the region o/t»(~ 70) in which the calculation is being performed. 

From the analysis in [13] we know that at large u ^ 1 the spectral function for lightlike momenta 
goes like ~ u'^^^ when m = 0. In [7], this behaviour seemed to persist for finite mass less than or 
equal to the critical one m < 1.3. These are the only stable black hole embeddings for vanishing 
D. To go to higher quark masses one must turn on D for stability. We have done a similar analysis 
and fitted our numerical results to a function of the form = a + f3vo'^ where the a, (3 and 7 are 
functions of m and D. The results are plotted in figure 5. The universality of 7 ~ 2/3 as observed 
in [7] corresponds to the the right most figure for values of m < rricrit ~ 1.3 where we see that 
all curves accumulate on an almost flat line. For masses approximately in the range m G (1.3, 6) 
the parameters of the fit exhibit a strong dependence on the baryon number density. This is to 
be expected, as these embeddings are the ones for which the induced horizon area changes most 
dramatically when D is switched on (see figure 2). For large m despite the fact that the embedding 
changes greatly in the UV, the change in the induced horizon area is negligible. For asymptotically 
large values of m the parameters seem to stabilize around some marginally D dependent values. 

4.1.3 Timelike momenta 

For timelike momenta both and X"*" are nonzero and contribute to the spectral function. In figure 
6 we have plotted both functions for several values of the baryon density, D = 0.1 and 0.5, in the 
range where the peaks are clear. As a function of the baryon density, the peaks seem to decrease in 
amplitude and increase in width monotonically. This corresponds to the poles moving further from 
the real axis in the complex plane. Concerning the spread among both components, the situation 
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m=4.54, q=20 

2500 V 




Figure 6: This plot shows the spectral function for longitudinal and transverse modes x"(c<-') and X'^i^) 
at q = 20 for two different values of D = 0.1, and 0.5. All curves oscillate around the zero temperature 
value given in eq. (1-3). 



is not easily captured from the figure, therefore we have tracked the movement of those peaks in 
figure 7. 



m=4.54, q=20 




Figure 7: Left hand plot: position of the first peaks in the longitudinal component x ■ Right hand side: 
difference between the positions of the v}^ peaks o/x" and x"*". All for different values of D = 0.1,0.25,0.5 
and 0.75. 

On the left we plot the position of the peaks in x", and see a drift towards higher values for 
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increasing D. On the right, we observe how the difference in peaks among both components indeed 
builds up also with D. 

In figure 8 we plot the difference in the full spectral function from the T = result as a function of 
to for q = 1 and different values of the baryon density. For values of m below (above) 1.3 the effect 
of the chemical potential is to enhance (suppress) the height of the oscillations as a function of t». 



m=l, q=l 



m=4, q=l 
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Figure 8: Finite T contribution to the full spectral function Ax{vo,q = 1) for different values of D. 



m=4 




5 10 15 20 



Figure 9: q dependence of the longitudinal spectral function at varying D and vo. 
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m=4, D=0.5 m=4, D=2 




Figure 10: For the transverse spectral function, the first four peaks move in the vo, q plane along 
hyperbolae with different limiting velocities, data points are marked with full circles. For the longitudinal 
modes we do not have a nice parametrization, so we plot only the raw data, marked as crosses. The thick 
line is the light cone q = to. 



Concerning the q dependence, figure 9 shows the longitudinal spectral function, x", for two different 
values of to = 10, and 20. For such high values we expect results very close to the zero temperature 
prediction as given in eq. (1.3) whose plot corresponds to the dashed line in the figure. 

We see peaks appearing in the region close to the light cone (q ~ to). As we move away from 
the light cone (to smaller q) the peaks in the spectral function for different tx> are in one-to-one 
correspondence. This could be seen by plotting the spectral function in the full t» — q plane. As 
before, we must treat peaks in the longitudinal and transverse spectral functions separately. We 
have plotted them in figure 10. Peaks in fit very well with a mass hyperbola of the form 

vo'~viy = M'^ (4.1) 
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where M„ should correspond to the modes of the quasiparticle spectrum, and presumably coincides 
with the continuation of the mass spectrum to the region where they become unstable (see the 
discussion in references [6, 44]). A careful analysis of the hmiting velocity, w„, of these unstable 
mesons was initiated in [45] and recently pursued in [33]. For D = 0.5, 2, 10 our data are consistent 
with the fits in table {A.2f. In [33] it was argued that the limiting velocity should be in any 
case v±^n = 1 for q — > oo. The interpolation to this behavior led the authors to speculate on the 
existence of an intermediate region where the group velocity may become superluminal (see [46, 47] 
for discussions on topics related to superluminal propagation). 





D 


= 0.5 


D = 2 


D 


= 10 


n 


M„ 




M„ 








1 


4.0 


1.0006 


4.2 


1.0038 


5.1 


1.01251 


2 


7.1 


1.0036 


7.8 


1.01086 


9.9 


1.02618 


3 


10.3 


1.0064 


11.5 


1.01762 


14.8 


1.03945 


4 


13.6 


1.0101 


15.2 


1.02422 


19.3 


1.04574 



(4.2) 



The departure from „ = 1 is certainly very small. However it displays a well defined monotonically 
increasing behaviour with D and m. Due to rotational invariance for q = we know that the curves 
for the transverse and longitudinal peaks must intersect at zero momentum. However, unlike the 
transverse case, we were unable to fit the longitudinal peaks to a mass hyperbola. 

4.2 Scalar modes 



m=l, D=2 m=4, D=2 




Figure 11: For larger values of the quark mass m the curves for the longitudinal spectral function and 
the scalar are more similar than for smaller values. 



^Wc fit the mass hyperbola eq. (4.1) by setting the value of Af„ equal to the peak position at q = and performing 
a fit analysis to the rest of the points. 
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The spectral function for the scalar modes is obtained along the same lines as stated in the previous 
section from the Gf^ component in eq. (3.9). Again, they must be calculated by solving the full 
coupled equation, although the operator becomes decoupled in the UV. In figure 1 1 we compare the 
longitudinal part of the spectral function and the scalar perturbation for two values of the quark 
mass m at fixed D. For large masses we recover the T = result and therefore the peaks narrow 
and become degenerate for the different modes, while for low quark masses we see that the peaks 
broaden, as expected, and move away from the longitudinal modes to which they are coupled in the 
JR. We see that the masses of the scalar quasinormal modes (the peaks in the spectral function) are 
higher than the equivalent longitudinal modes. The lifetimes are also shorter, given by a widening 
of the scalar peaks. Note that this same pattern is seen in a slightly different context in [48] for 
instance. 



5. Conductivity 

The electrical conductivity can be obtained from the zero-frequency slope of the spectral function 
for electromagnetic currents. In [7] the conductivity of the plasma was calculated. We may do the 
same here, however we include the effects of a finite baryon density. This quantity may be calculated 
at null-momentum [13, 7] and therefore only the transverse spectral function contributes; 



(5.1) 

fcO = |k| 



We show the results for this numerical calculation in figure 12. In [23] the conductivity of the 
plasma was obtained from Ohm's law using AdS/CFT techniques. In the presence of an electric 
field e a current establishes satisfying Ohm's law J = ae. The electrical conductivity is given by 
the following expression ^: 

^=^^v/^^(l-^g)=^ + ^, (5.2) 

where e the electric field strength, is zero in the present setup. An insightful interpretation of 
this formula was given in the aforementioned reference in terms of two sources of charge carriers 
contributing to the conductivity. One comes explicitly from the baryon charge asymmetry, and is 
given by the contribution. The other comes from charge carriers produced via pair creation in 
the plasma. This effect clearly depends upon the mass of the fundamental charged quarks through 
the combination (1 — -^q) and is equal to 1 for massless quarks and to for infinitely massive ones. 
This second effect is implicitly dependent on D. Our numerical results plotted in figure 12 match 
perfectly with the above analytic expression in the limit of vanishing electric field 5 — 0. 
This result can also be derived from the exact analytic expression for the low frequency limit of 
the spectral function as given by (C.6). The calculation in [23] is a macroscopic result. The fact 
that this matches with our derivation from microscopic linear response theory is a highly non-trivial 
check of consistency since both explicit and implicit dependence on the baryon density are correctly 
reproduced from the full spectral function. 

■^Note that our definition of a and D differs from that in [23] by constant coefficients which are given in detail in 
appendix E. 
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Conductivity 



D=4 

1.0 - 



0.8 - 
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D=2 



0.4 - 




Figure 12: Conductivity, calculated numerically from the derivative of the spectral function as a function 
of the quark mass for various values of the baryon number density. The plot matches with the analytic 
expression given in eq, (5.2) upon setting e = 0. 



6. Conclusions 

Exploiting the AdS / CFT correspondence at finite temperature and baryon density we have consid- 
ered a holographic dual of = 4 super Yang Mills with quenched massive flavours (which break 
supersymmetry to A/" = 2). These flavours are constructed from a small number, Nf, of D7-branes 
embedded in the black hole geometry sourced by a stack of D3-branes. The addition of finite 
baryon number is modeled by turning on a vev for the worldvolume gauge field Aq living on the 
D7-brane. Working at fixed baryon number gives an additional parameter D in addition to the con- 
stituent quark mass in units of the Hawking temperature. This opens up the possibility of studying 
the spectral functions for flavour bilinear operators in regimes where there is clearly recognizable 
structure. This structure is interpreted in terms of quasinormal modes of mesons which develop 
a finite lifetime, and therefore a broadening of their spectral functions at finite temperature. By 
looking at the spectral function of longitudinal and transverse modes as well as excitations of the 
scalar embedding coordinates on the D7-brane we have studied the behaviour of the poles in the 
spectral function as a function of D as well as spatial momentum. In particular, for finite D and 
q we are forced to solve for the scalar and longitudinal modes as a coupled system of differential 
equations. 

In general we find a splitting of the peaks between the spectral functions X"*" and x" which increases 
with both the baryon number and the spatial momentum. Plotting this effect for different quark 
masses, we find that the effect is most pronounced for values of m in the range ~ (1.3,4). This 
is where the induced horizon area changes most upon switching on D. We have also investigated 
the dispersion relations i»(q), where t» is the position of the peaks in the spectral function. For the 
peaks in x'^ "we find a very good fit of the data for the position of the peaks on a mass hyperbola 
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of the form to = ^/v'^<\'^ + M^, as discussed in [33]. For low masses, close to the critical mass, (as 
defined by the phase transition present in the absence of baryon density) and very low values of D 
(so that the peaks are still clearly seen) the values of v±^n ^ 1 found in [33] are reproduced. For high 
enough values of the baryon density (and therefore high quark masses, in order to keep the position 
of the poles close to the real axis) we also find values of v±^n > 1- For the longitudinal modes we do 
not know the correct function to fit. All of these issues should be confronted with the observation 
made in [33] that as q ^ oo the zero temperature relativistic dispersion relation with t» = q should 
be recovered. The clarification of this issue demands the use of powerful techniques, in particular 
the analysis of the quasinormal modes. Certainly this will be a very challenging calculation for the 
coupled longitudinal and scalar modes, where the technique of converting the equation of motion 
to Schrodinger form becomes very difficult. We will come back to this issue in the future. 
The most spectacular effects of the baryon density appear in the small xx> region for lightlike 
momenta. This is seen in the set of figures 3 where most of the change with D occurs for to < 2. 
There are two main consequences of this. On one hand the baryon density is responsible for an 
increase in the photoproduction rate, for small quark masses by a very large percentage. On the 
other, the slope as a function of a; at the origin is directly related to the electrical conductivity. 
We have provided analytical and numerical evidence that the effect of baryon number is given by 

a ~ \J {1 — iPq)^ + D'^. This provides a highly nontrivial check of the results obtained in [23] for 
the electrical conductivity in a macroscopic setup. 

It is clear that there is room for more research in this area, in particular a study of the limiting 
velocity and functional form of the dispersion relations for the longitudinal and scalar modes would 
be of great interest. In order to investigate this accurately, a large amount of computational time 
will be necessary. The study of the effects detailed in this paper may also be applied to various other 
finite temperature geometries. As mentioned above, a full calculation of the quasinormal modes 
will allow us to go to higher values of spatial momentum than are currently available by studying 
peak positions in the spectral functions. 
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A. Equations of motion 

In this section we shall give more details about the equations we have been dealing with in the text. 
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We will try to keep quite general and provide formulae that encompass general Dp/ Dq intersections. 
The background metric will be given by a type II supergravity solution, of the following generic 
form 

ds]^ = GMN{X)dX^^dX^ , M,iV = 0,...,9. (A.l) 

We shall consider the embedding in the following gauge: separate into two groups X"", a = 
0,l,...,g and I = g + 1,...,9. 

fX'^ = e = {xO,x%«,^n (a = 0,l,...,g) 1 

\x^ = Z^(0 = M^n (/ = g + l,...,9) J ■ ^ 

Then the pullback metric and the dilaton are 10 dimensional bulk fields in origin. They are evaluated 
on the Dq brane 

= 0(^,^(0) ; gab = G,,{^,Z{0) + Gjj{^,Z{0)^^, (A.3) 

where we have assumed that the background metric has no mixed terms Gai = 0. On the other 
hand, the gauge fields are purely world-volume fields, hence 

Fa,{0 = daAtiO - dMO ■ (A.4) 
Let us consider a set of perturbations of the form 

z^iO = + ezW'iO ; MO = + , (A.5) 

and expand up to second order in e. Naming gab the pullback metric we have 

= 0W + 60« + 6 V^) ; 9ab = 9iS + eg^ + e^g^ ; = + eF^ . (A.6) 

For example 

0(1) = 9,0(°)z(i)^ ; 0(2) = ^didj<P^'h^''>-'z^'^' , 

nW-fr r(o)/^(o)^\ ^(i)K , 9^ Jo)i (i)j 

A'^)_^(r 4-r J{^)\ Jl)K Jl)L , I Jl)K(0)I {l)J 

9ab — 2 \'-'ab,KL + <^IJ,KLZla ^,b J ^ ^ +'-^IJ'^,a ^,b +^'-^IJ,KZ' Z(„ Z^^s^ , 

We want to expand the DBI lagrangian 



C = e-W\detig + F)\, (A.7) 
in powers of e. Note that the absence of the Wess-Zumino term is detailed in appendix D. Defining 

_ -^(0) _ (0) p(0) . Ml) _ , p(l) . ^(2) _ ^(2) . . 

lab = yab - 9ab + ^ ah ' ^ab " 9ab + ^ ab ' ^ab " 9ab ' l^'^J 
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as well as 

the lagrangian up to C(e^) is 



cb 1 



C = e- VI det g\=Co + tCi + e^A , (A.9) 

with 7 = — det we have 

A = e-^<"'v^(-0(^) + ltrS«)) , 

£2 = e-*^"'v^ f JtrE(2) - itrE(^)2 + ^(trS^i))^ - l0(i)trS(i) + U^'^' - . (A.IO) 

y z 4 o 2 2 y 

Let us now focus on a particular family of bulk metrics adapted to the setting of a Dp/Dq inter- 
section, with q = p + n + 1. 

dsji^B = Goo{u)dxl + Gii{u)dXp + Guu{u)du^ + Gee{u, 'ip)d^ll + G^^{u, %l))d%l)^ + G^^{u, ilj)dQ.'^_p_^ , 

Here cifi^ is the metric on a unit n dimensional sphere wrapped by the flavour brane. Similarly the 
perpendicular space has been written in adapted polar coordinates {tp,ip^). We will be considering 
perturbations of the following form 

AaiO ^ SaoAoiu) + ee-^('^-'"-^^") , (A.ll) 

Notice that with this ansatz, we have (p^'^^u) and (j)^^'' = (p^"^^ = in (A. 6). The equations of 
motion for the background profiles Aa{u) are obtained from Cq and give 

^ ^ 2du (e-^'°' V7^'7""G'^V') - e'^"" V7 (^"7""^^^,^ + n^'Gee,^) = , (A. 12) 
^0 ^ 5.(e-'^'"V77°") =0. (A. 13) 

From C2 we obtain the equations of motion for the perturbations. The transverse fluctuations 
A±{u) = .42,3 (m) satisfy a decoupled equation 

EoM[A^] ^ A'[ + du log [e-f'"' v^7''7""J ^1 - ^ ^ ^± = . (A.14) 
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On the other hand the longitudinal fluctuations couple to the proflle perturbation 



EoM[Ao] ^ X + d^\og e-^^"^7°V" X - q^{qAo + cuA,) + q'^^J'G^^^ 



11 



-.2 1 / „iJ, 



7 



Ou 



I I 



(A.15) 



^00 



(A.16) 



EoUm ^ *" + a„log e-'^^°'v^7""G'v,v,(l-V''A) + 



A(l - V^'A) 



7«"(1 - V''A) 

and there is one constraint from the gauge choice, Au — 



^/ 11 On 7/;S°" 



Ou 



where 



(A.18) 
(A.19) 



and ii{u) is a rather lengthy expression we give here for completeness 



H{u) 



du 



7""(1 - V''^7""G'VV') 



2r^G^fil - ^'27""G^v) 
_ (V^-^g^^,,,,-|7""(g^^,^)'V^^^) .J. 

We work with the gauge invariant field combination Z = ~i e^i'^^o-qxi)pW _ -)-^^^. Prom the 
equations for these gauge perturbations we see that taking combinations 

q EoM[Ao] + u EoM[Ai] + iquj \[ ^ EoM[Au] , 



log' (e-'^""v^7"7"") - 



A'-S' 



(A.21) 

EoM[Au] , 
(A.22) 
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allows us to write the equations of motion in terms of Z and ^ as follows 



Z" + %^Z' + + Ci*" + + = 
Z" + %2Z' + :s&2Z + + 5?2*' + = 



(A.23) 
(A.24) 



with 

Iti = log' 



112 = log' e-^'°'V7yi7"- 

•Si = -^2 



2^,11 



^2^00 + ^2^11 



A' - S 



^2^00 + ^2^11 



„0m 



£1 = -X!—^ 



55i = -4?^(S + A')- 



■''"a log' 



^2 



log' e-<^™V77""G'^^(l-V''A) 



00 



g7 



On 



(^2^00 _|_ g2^,ll ^^,,00^^,UM 



log' (e-^'°'V77^^7''") 



97 



Ou 



;' / 



«2 



97 



^00 ^uu 



S'- A 

00 



^2 ^00 ^ ^2^11 



go; 



^2^00 + ^2 



A'-S 



log 



^.2 7-^ 

7"" ^2^00 _|_ ^2^11 



5 (log' (e-^"'V77V'^) - 



7~ 
A'-S 



g(l-^'A) 



In the particular case of the D3/D7 intersection, we can express the factors governing the equations 
of motion as 

7oo = -/H — z. — ; 711 = 722 = 733 = 



u 



u 



In 



L2(1 - V2 + 4mV(m)^ 
4(1 - ^jj^)u^f{u) 



; 7o« = -7uo = -27ra'/lo(M) , 
L2 



(A.25) 



1-^2 



n = 3 



0. 
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with inverse components 



7°° = ; = 7^^ = 7=*^ = - ; r = -. (A.26) 

7oo7«M + (7o«) 7ii lee 

^uu _ Too . ^Ou _ _^uO _ ~70n 

7oo7w« + (7on)^ ' 7oo7«« + (7o«)^ ' 

B. Boundary conditions for the coupled system 

In this paper we deal with a coupled system of second order differential equations with a regular 
singular point at p = 0. That is 

z"(p) + ^z-(p) + ^z(p) + ^np) + ^m = o, (B.i) 
*"(P) + ^np) + + ^z'(P) + ^z(P) = , (B.2) 

p p^ p p^ 

where the singularity at p = is such that all the functions {A{p), B{p), . . . C{p), D{p)} are regular 
and can be Taylor-expanded in powers of p. This expansion will be indicated in general as 

oo 

M{p) = Y,m^P\ (B.3) 
i.e., using small letters. Performing a Frobenius expansion for Z and \1/ with indices A and r] 

oo oo 

z = p'Y,^^p' ; ^ = p"5Zv^^pS (b.4) 

a system of recursion relations for the coefficients Zi and ipi can be calculated 

oo k 



\J2J2\P^ + i)(A + i - l)6ik + (A + i)ak-i + bk^i) Zi + p" {{r] + i)ck-i + dt-i) p'^ = 

fc=0 1=0 
^ oo fc 

- ^ ^ [p'' [ir] + i){T] + i- l)6ik + {r] + i)dk-i + h-^ + ((A + i)ck-i + 4-i) z. 



We can now focus on the D3/D7 background. Close to the horizon we take p = 1 — u. For equations 
(A. 24) and (A. 23) we obtain the following coefficients 

ao = fio = 1 ; bo = bo= ^^^^^^ ; cq = cq = ; do = do = , (B.5) 
and find two solutions of the indicial equation 

Ai = = V± . (B.6) 
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There is one free parameter, say zq (or ipo). Now for an analysis close to the boundary let p = u in 
which case, we get 

ao = , bo= co = do = 

; 7 3 ; . (B.7) 

ao = -1 ^0 = ^ co = do = 0, 

giving the following exponents 

Ai,2 = {0, 1} ; ^1-2 ={^,^}. (B.8) 

With these coefficients we find near the boundary 

Z{u) =A Z^^^ {u)+uB Z(^) (u) , (B.9) 

^{u) =u'/^M¥^'\u)+u^^^C¥^\u), (B.IO) 

where 

oo 2 2 ^ 

Z(^\u) = 1 + ^24"^^' -^^^^og uZ(^\u) ■ Z(^)H = 1 + X:4V, (B.ll) 

i=l i=l 

OO 2 2 ^ 

^W(m) = i + yz^;-V-^^i^MiogM^(^)(M) ; ¥''\u) = i + y4'\\ (b.i2) 

i=l 1=1 

C. Asymptotic expressions for lightlike momenta 

In order to trust the numerical results it is important to compare both to previous results, or 
to analytical expressions wherever possible. In the current setup it is possible to calculate such 
analytic results only in certain asymptotic limits. Setting the momenta on the light cone, t» = q, 
the longitudinal polarization vanishes and the equation for the transverse component of the gauge 
field (A. 14) acquires the following form upon introducing the particular values given in (A. 25) and 
(A.26) 

Al^g^,^(f:p.^I^] ^^^^n^zIll^A.^,. (C.I) 

In contrast to the massless, D = {) case [13], we have found no analytic solution when baryon number 
is turned on. We can however perform analytic expansions in the small and large xo regimes and 
extract analytic information perturbatively. 

C.l Low frequency limit 

For the small xo limit we will use a perturbative approach to find the analytic behaviour. We begin 
extracting the regular singularity at m = 1 by substituting 

A^{u) = r^/' (Yoiu) + tt)Fi(n) + 0{w')) , (C.2) 
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with Yn{u) regular at m = 1. We can then expand the equation of motion as a series in to and solve 
it order by order. From (C.l) we are led to the following equation for Yo[u), 



fip\l ip(^ + D^u^ 

Clearly to avoid a singular solution at m = 1 we must set Ci = and hence Yq is a constant. For 
Yi we perform the same procedure, finding the general solution 



which is regularized at the horizon setting 

C2 = iYo^{l-ijir + D'^, (C.5) 

where ipo = = !)• 

We can calculate the spectral function by looking at the m — > end of the solution and find that 

NrNfT^ 



X", - ^-^^ro^il - + D' + 0{w') , (C.6) 

which leads to precisely the same result for the conductivity as that obtained in the vanishing 
electric field limit in [23] . This is a very elegant result showing that the microscopic and macroscopic 
calculations agree with each other. 

C.2 High frequency limit 

Now we move to the t» ^ 1 limit. We use the Langer-Olver method (a version of the WKB approx- 
imation) to construct the asymptotic solution and will consider only the massless case. Following 
[13, 20], we perform the following transformation 



^±(^) = \/y^ ^^;p2^3 y{u); u = -x, (C.7) 

for the equation (C.l), in order to rewrite it in a Schrodinger form 

y"{x) = [v0'H{x) + GixMx), (C.8) 

where 

^' xf(xy- l-DV ' ^ ' 

X = —u G [—1,0], and G{x) is another function which will not be needed in the following. For 
large to the dominant term on the right-hand side of (C.8) has a simple zero at x = and thus. 
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according to [49] , the asymptotic solution can be expressed in terms of Airy functions. With this 
goal in mind, we change variables from x to ( as 



c 



dC 

dx 



H(x) 



c 



2/3 



(C.IO) 



and rescale y{x) to W{Q as 



y 



dC 
dx 



-1/2 



w. 



After which we are left with a new differential equation for W{Q 

d^W 



(C.ll) 



(C.12) 



where 77(C) is another function whose form is not needed in the analysis^. For large t» (C.12) reduces 
to Airy's equation and the solution to leading order is 



W{C) = AoAi(w2/3() + 5oBi(t»2/3() + . 



(C.13) 



Dictated by the incoming-wave boundary conditions at the horizon we set Bq to zero. Thus the 
approximate solution for A±{u) for large t» is 



AAu) 



An 



a-u) 



1/4 



leading to the following expression for the transverse correlator 



NcNfT^ 



Pi 
5 



^_l)2/33l/3r(2/3) 



2/3 



(C.14) 



(C.15) 



Note that the transverse correlator depends on D, but only through the real part. The trace of the 
spectral function in the high-frequency limit for lightlike momenta is ^ 



NcNfT^ 3^/6r(2/3) 
2 21(1/3) 



to 



2/3 



0(t») 



(C.16) 



and is independent of the presence of baryon density, at least to leading order in to, giving the same 
result as in [13]. The appearance of the baryon density in the real part of the transverse scalar 
indicates possible deviation of our result from that of [13] at the first subleading order in to. The 
numerical solution matches this claim. 

■*The exact expression for this function in terms of G{x) and H{x) is in [20] 

^In the niassless case the differential equation for the scalar mode is decoupled from the longitudinal one, and 
its contribution to the spectral function can be calculated independently. In the lightlike case, even if the initial 
differential equation for the scalar is different from that for A±_ , we can prove that the contributions to leading order 
in Lo in the transverse scalar and consequently to the spectral function are (C.15) and (C.16). 
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D. The Wess-Zumino term 



In order to illustrate that there is no contribution from the WZ term to the equations of motion, 
we start by writing the Ramond-Ramond four-form C4 in the following notation: 



-dt A dxi A dx2 A dx^ . 



The contribution of this four-form to the five-form field strength is: 



dCA 



-dt A dxi A dx2 A dx^ A du . 



whereas the Hodge-dual contribution is: 

i^dCi = 4L^(1 - 9 cos 9 dip A d(j) A d9 A d02 A d(j)s . 

From this expression we can deduce the four-form dual to the one in (D.l): 

2L^ 



(D.l) 

(D.2) 
(D.3) 



+ 



3 



- (sin^ 9dilj Ac/0 - sin 29 (p dip A d9) 

iP^' 



sin 29d(p A d9 



A 



(D.4) 



P 



Next, we use (D.4) to calculate its pullback onto the D7-brane, so we have: 

= —{l-ip^)ip sin^ 9p^^dx'' A dx^ A d<P2 A 
3 ^ ox"' oxP 

2L^ dib 

+ {l-ip'^)ip<p sin 29-— dx"" Ad9 A d<p2 A d<p^ 

3 ox°' 



Ip"' 



sm29-—dx'' Ad9 A 



A d(p3 . 



(D.5) 



2 J 

We are interested in WZ terms with a product of two fields which are fluctuations on the D7-brane. 
Each term in the above four-form will contain at least one factor of a fluctuating field. The pull- 
back must be contracted with two factors of the brane world-volume gauge field, and therefore one 
of these two factors must be the background component Fqu corresponding to the presence of the 
finite baryon number density. The other component will therefore be a fluctuation of a gauge field. 
The eight-form product which will contribute at second order will have the following three terms: 



F AFAP 



Fou AFuAi L'lP^ { 1 



dx^ Ad9 A 



A 



+ -fou A -^23 

+ Fo„ A F31 (l'^iP^ ( 1 



y) sm29^^]dx' Ad9A 



sm29^]dx' Ad9A 
2 / ox^ 



A 



A 



(D.6) 
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This can be written as: 

^ r d(j) 1 

1A2 - d2Ai)- h cyclic in 1,2,3 Yoim ■ 

0x3 J 

Then it is easy to produce the WZ contribution to the equation of motion for 0: 



-A'qL^^^ ( 1 - y j sin 2^ 



(D.7) 



-A'^L'^ (^1 - "ipj sin 26 [ds {diA2 - ^2^1) + (^2^3 - ^3^2) + 82 {d^A, - d^As)] . (D.8) 

From this equation it is obvious that the terms coming from the WZ part of the action vanish due 
to the antisymmetry of the gauge field indices. 

For the field Ai the WZ contribution to the equation of motion is given by: 



_d_ 



_d_ 



A'L'i, I 1 



sin 26 



dxi 



(D.9) 



which once again gives no contribution from the WZ term. For the fields A2 and ^3 the above 
argument goes through identically. 



E. Comparison with the conductivity found by Karch and O'Bannon 

In [23] the following result for the conductivity was obtained: 



\ 



^ C0S6^ ' 



The difference in normalisation between our conductivity, a, and theirs, a, is: 

4 NfN^T 



(E.l) 



4:71 



(E.2) 



The factor ^^^"'^ comes from the normalisation of x ^^id the factor 4 is a difference in definition 



between our conductivities. Therefore: 



(1 - ^lY + 



167r2 L)2 



Inserting their definition of A and remembering to put back in the AdS radius gives: 



(E.3) 



a= -^l{l-i,lY + 



(E.4) 
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The difference between our definitions of D is calculated by looking at the UV behaviour of Aq. We 
must also note that the difference in our definitions of the radial coordinate are given by: 



(E.5) 



where zh = 

This gives a difference in Ds as: 



NfNrT^L'^ 



where again the tilded quantity is their definition. Using this in the definition of the conductivity 
gives us: 

a=j^J{l-ror + D^. (E.7) 
This result matches perfectly with our numerical calculation given in section 5. 
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